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Abstract  Three  methods  have  been  used  to  experimen¬ 
tally  determine  the  roughness  function  (A17+)  for  several 
rough  surfaces.  These  include  the  rotating  disk,  the 
towed  plate,  and  the  velocity  profile  methods.  The  first 
two  are  indirect  methods  in  as  much  as  they  rely  on 
measurements  of  overall  torque  or  resistance  and 
boundary  layer  similarity  laws  to  obtain  AU+,  whereas 
the  velocity  profile  method  provides  a  direct  measure¬ 
ment  of  AfJ1".  The  present  results  indicate  good 
agreement  between  the  towed  plate  and  the  velocity 
profile  methods  for  all  of  the  surfaces  tested.  Tests  for 
the  rotating  disk  were  carried  out  at  much  higher  unit 
Reynolds  numbers.  Using  this  method,  the  results  for 
sandpaper  rough  surfaces  agree  within  their  uncertainty 
with  a  Nikuradse-type  roughness  function  in  the  fully 
rough  regime,  while  a  spray  painted  surface  agrees  with  a 
Colebrook-type  roughness  function. 

Nomenclature 

B  smooth  wall  log-law  intercept,  =5.0 

CF  overall  frictional  resistance  coefficient, 

=  (ft  )/(\pu}s) 

Cf  skin-friction  coefficient,  =  (t0)/ (|pU2) 

Cm  torque  coefficient,  =(2M)/(pR5(0eo)2) 

Fd  drag  force 
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arbitrary  measure  of  roughness  height 
acceleration  parameter,  = 
plate  length 
torque 

number  of  samples  or  replicates 

disk  radius  N 

centerline  average  roughness  height,  =  jjYl  bi'l 


root  mean  square  roughness  height,  =  y^Xjy2 

maximum  peak  to  through  height,  =ymax-ymin 

5 

ten  point  roughness  height,  =  X]  Omaxi  ~  y mini) 

i=i 

Reynolds  number  based  on  plate  length,  =UeL/v 

Reynolds  number  based  on  disk  radius,  =</>ojR2/v 

momentum  thickness  Reynolds  number,  =Ue0/v 

wetted  surface  area 

mean  velocity  in  the  x  direction 

freestream  velocity  relative  to  surface 

friction  velocity,  =  y/ z0/p 

streamwise  fluctuating  velocity 

wall-normal  fluctuating  velocity 

streamwise  distance  from  plate  leading  edge 

normal  distance  from  the  boundary 

boundary  layer  thickness 

roughness  function 

roughness  function  slope,  =d(AU+)/d(ln/c+) 
wall  datum  offset 

swiri  factor,  (j.)  J 
von  Karman  constant,  =0.41 
kinematic  viscosity  of  the  fluid 
wake  parameter  5  2 

momentum  thickness,  =  f  (l  —  ^f)  d y 
density  of  the  fluid  0 

wall  shear  stress 
angular  velocity 


Superscript 

+  inner  variable  (normalized  with  Uz  or  UJv) 


Subscripts 

en  disk  rotating  in  enclosed  tank 

min  minimum  value 

max  maximum  value 

R  rough  surface 

S  smooth  surface 

00  disk  rotating  in  unconfined  fluid 
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Introduction 

Wall-bounded  turbulent  flows  over  roughness  occur  in  a 
multitude  of  situations  in  both  the  natural  environment 
and  in  engineering  applications.  Examples  range  from 
wind  over  plant  canopies  to  boundary  layers  on  ships  and 
aircraft.  Since  surface  roughness  affects  not  only  the  near 
wall  mean  velocity  profile  but  also  the  resulting  drag  and 
heat  transfer,  it  has  been  the  subject  of  a  great  deal  of 
research  (e.g.,  Nikuradse  1933;  Clauser  1954;  Hama  1954; 
Ligrani  and  Moffat  1986;  Raupach  et  al.  1991;  Antonia  and 
Krogstad  2001). 

The  mean  velocity  profile  in  the  inner  region  of  the 
smooth  wall  turbulent  boundary  layer  is  given  by  the 
classical  log-law. 

U+  =  -ln(y+)  +B  (1) 

K 

Clauser  (1954)  demonstrated  that  the  primary  effect  of 
surface  roughness  was  to  cause  a  downward  shift  in  the 
log-law.  This  downward  shift,  ALT'",  is  called  the  roughness 
function.  The  generalized  log-law  for  rough  and  smooth 
walls  is  therefore  given  as 

U+  =  —In  (y  +  e)+  +  B  —  AU+  (2) 

K 

For  so-called  “k-type”  roughness,  A l/"1"  is  a  function  of 
the  roughness  Reynolds  number,  k+,  defined  as  the  ratio  of 
the  roughness  length  scale,  k ,  to  the  viscous  length  scale, 
v/UT.  Once  AU+=AU+(k+)  for  a  given  rough  surface  is 
known,  it  can  be  used  in  a  boundary  layer  code  or  a 
similarity  law  analysis  to  predict  the  drag  of  any  body 
covered  with  that  roughness  (Granville  1958).  However, 
determining  the  roughness  function  is  not  without  diffi¬ 
culty.  There  are  two  primary  reasons  for  this.  First,  the 
roughness  function  is  not  universal  for  all  roughness 
types.  Nikuradse’s  (1933)  pipe  flow  experiments  with 
closely  packed,  uniform  sand  indicate  this  roughness  type 
has  a  universal  roughness  function  with  k  simply  being  the 
diameter  of  the  individual  sand  grains.  However,  the 
roughness  functions  for  most  naturally  occurring  surfaces 
do  not  behave  like  closely  packed,  uniform  sand.  Cole- 
brook  (1939)  demonstrated  this  in  a  study  of  the  irregular 
surface  roughness  in  pipes  resulting  from  the  manufac¬ 
turing  process.  The  other  primary  difficulty  is  specifying  a 
proper  roughness  length  scale,  k,  for  a  given  roughness 
geometry.  Granville  (1987)  states  that  more  than  one 
roughness  length  scale  may  be  required  to  adequately 
represent  an  arbitrarily  rough  surface.  For  these  reasons, 
determining  the  roughness  function  for  a  generic  surface 
requires  experimental  tests  (Grigson  1992). 

Various  experimental  methods  have  been  offered  for 
determining  AU1"  for  rough  surfaces.  The  most  evident  is 
to  directly  measure  A  U+  using  the  mean  velocity  profile  in 
the  log-law  region  over  the  roughness  (Clauser  1954). 
However,  the  determination  of  UT  for  a  rough  wall  profile 
is  more  prone  to  error  due  to  the  two  additional  param¬ 
eters  (the  position  of  they  origin  and  A  U+)  that  must  also 
be  found  when  fitting  the  profile  (Acharya  et  al.  1986). 
There  are  several  indirect  methods  for  finding  ALU  as  well. 


These  include  measuring  the  pressure  drop  in  fully 
developed  pipe  or  channel  flow  (Nikuradse  1933),  the 
torque  on  rotating  disks  (Granville  1982),  and  the  drag  on 
towed  flat  plates  (Granville  1987).  Using  these  measure¬ 
ments  and  boundary  layer  similarity  laws,  ALT'"  can  be 
determined.  The  goal  of  the  present  experimental  inves¬ 
tigation  is  to  compare  the  roughness  functions  obtained 
using  three  of  these  methods;  the  velocity  profile,  rotating 
disk,  and  towed  plate  methods. 

2 

Experimental  methods  and  conditions 

In  this  study,  four  surfaces  were  tested  for  each  of  the 
roughness  function  determination  methods.  These 
included  a  baseline  smooth  surface,  a  60-grit  sandpaper 
surface,  a  220-grit  sandpaper  surface,  and  a  sprayed  epoxy 
surface.  The  surface  profiles  of  the  rough  surfaces  were 
measured  using  a  laser  diode  point  range  sensor  (model 
No.  PRS  40,  CyberOptics  Corp.,  Minneapolis,  Minn.,  USA) 
laser  profilometer  system  mounted  to  a  Parker  Daedal 
(model  No.  106012BTEP-D3L2C4M1E1,  Parker  Hannifin 
Corp.,  Irwin,  Penn.,  USA)  two-axis  traverse  with  a  reso¬ 
lution  of  5  pm.  The  resolution  of  the  sensor  is  1  pm  with  a 
laser  spot  diameter  of  10  pm.  Data  were  taken  over  a 
sampling  length  of  50  mm  and  were  digitized  at  a  sam¬ 
pling  interval  of  25  pm.  Ten  linear  profiles  were  taken  on 
each  of  the  test  surfaces.  The  sandpaper  surfaces  were 
constructed  from  the  same  sheet  of  sandpaper,  so  the 
roughness  on  these  surfaces  agreed  within  their  uncer¬ 
tainty.  The  sprayed  epoxy  surface  used  in  each  method 
varied  somewhat  because  of  differences  in  paint  applica¬ 
tion.  The  mean  roughness  statistics  for  each  of  the  test 
surfaces  along  with  their  95%  precision  confidence  limits 
are  given  in  Table  1. 

2.1 

Velocity  profile  method 

The  velocity  profile  method  experiments  were  carried  out 
in  the  closed  circuit  water  tunnel  facility  at  the  United 
States  Naval  Academy  Hydromechanics  Laboratory.  The 
test  section  is  40  cm  by  40  cm  in  cross-section  and  is 
1.8  m  in  length,  with  a  tunnel  velocity  range  of  0-6.0  m/s. 
In  the  present  study,  the  freestream  velocity  was  varied 
between  ~1.0  m/s  and  3.5  m/s  (Rex=1.4xl06-4.9xl06). 
Flow  management  devices  include  turning  vanes  placed  in 
the  tunnel  corners  and  a  honeycomb  flow  straightener  in 
the  settling  chamber.  The  honeycomb  has  19  mm  cells  that 
are  150  mm  in  length.  The  area  ratio  between  the  settling 
chamber  and  the  test  section  is  20:1,  and  the  resulting 


Table  1.  Roughness  statistics 


Specimen 

Smooth 

Ra 

(pm) 

NA 

Rq 

(pm) 

NA 

Rt 

(pm) 

NA 

Rz 

(pm) 

NA 

60-Grit  sandpaper 

126±5 

160±7 

983±89 

921±82 

220-Grit  sandpaper 

30±2 

38±2 

275±17 

251 ±14 

Epoxy — velocity  profile 

9±1 

12±1 

76+8 

71±7 

Epoxy — rotating  disk 

1 5±  1 

19+1 

128±13 

120± 13 

Epoxy — towed  plate 

10±1 

17±1 

135±7 

130±6 

Uncertainties  represent  the  95%  confidence  precision  bounds 
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Fig.  1.  Schematic  of  the  flat  plate  test 
fixture  for  the  velocity  profile  method 


freestream  turbulence  intensity  in  the  test  section  is 
-0.5%. 

The  test  specimens  were  inserted  into  a  flat  plate  test 
fixture  mounted  horizontally  in  the  tunnel.  The  test  fix¬ 
ture,  shown  in  Fig.  1,  is  similar  to  that  used  by  Schultz  and 
Flack  (2003).  The  fixture  is  0.40  m  in  width,  1.68  m  in 
length,  and  25  mm  thick.  It  is  constructed  of  a  high  den¬ 
sity  foam  core  covered  with  carbon  fiber  reinforced  plastic 
skins  and  was  mounted  horizontally  in  the  tunnel’s  test 
section  along  its  centerline.  The  leading  edge  of  the  test 
fixture  is  elliptically  shaped  with  an  8:1  ratio  of  the  major 
and  minor  axes.  The  forward  most  200  mm  of  the  plate  is 
covered  with  36-grit  sandpaper  to  trip  the  developing 
boundary  layer.  The  use  of  a  strip  of  roughness  was  shown 
by  Klebanoff  and  Diehl  (1951)  to  provide  effective 
boundary  layer  thickening  and  a  fairly  rapid  return  to  self¬ 
similarity.  The  test  specimen  mounts  flush  into  the  test 
fixture  and  its  forward  edge  is  located  immediately 
downstream  of  the  trip.  The  removable  test  specimens  are 
fabricated  from  12  mm  thick  cast  acrylic  sheet  350  mm  in 
width  and  1.32  m  in  length.  The  boundary  layer  profiles 
presented  here  were  taken  1.35  m  downstream  of  the 
leading  edge  of  the  test  fixture.  Profiles  taken  from  0.75  m 
to  the  measurement  location  confirmed  that  the  flow  had 
reached  self-similarity.  The  trailing  150  mm  of  the  flat 
plate  fixture  is  a  movable  tail  flap.  This  was  set  with  the 
trailing  edge  up  at  -5°  in  the  present  experiments  to 
prevent  separation  at  the  leading  edge  of  the  plate.  The 
physical  growth  of  the  boundary  layer  and  the  inclined  tail 
flap  created  a  slightly  favorable  pressure  gradient  at  the 
measurement  location.  The  acceleration  parameter  ( K ) 
varied  from  7.4x1 0-8  at  the  lowest  freestream  velocity  to 
2.0xl0~8  at  the  highest  freestream  velocity.  The  pressure 
gradient  did  not  vary  significantly  between  the  test 
specimens. 

Velocity  measurements  were  made  using  a  TSI  IFA550 
(TSI  Inc.,  St.  Paul,  Minn.,  USA)  two-component,  fiber¬ 
optic  LDV  system.  The  LDV  used  a  four-beam  arrange¬ 
ment  and  was  operated  in  backscatter  mode.  The  probe 
volume  diameter  was  -90  pm,  and  its  length  was 
-1.3  mm.  The  viscous  length  (v/Uj  varied  from  a  mini¬ 
mum  of  5  pm  for  60-grit  sandpaper  at  the  highest 
Reynolds  number  to  24  pm  for  the  smooth  wall  at  the 
lowest  Reynolds  number.  The  diameter  of  the  probe  vol¬ 
ume,  therefore,  ranged  from  3.8  to  18  viscous  lengths  in 
the  present  study.  The  LDV  probe  was  mounted  on  a 
three-axis  traverse  unit  (Velmex  Inc.,  Bloomfield,  N.Y., 
USA).  The  traverse  allowed  the  position  of  the  probe  to  be 
maintained  to  ±10  pm  in  all  directions.  In  order  to 


facilitate  two-component,  near  wall  measurements,  the 
probe  was  tilted  downwards  at  an  angle  of  4°  to  the  hor¬ 
izontal  and  was  rotated  45°  about  its  axis.  The  tilting  of  the 
LDV  probe  led  to  a  small  bias  error  in  the  V  measurements 
because  of  the  introduction  of  the  w  fluctuations  in  the 
measurement  plane.  The  error  in  the  measured  Reynolds 
shear  stress  caused  by  tilting  was  <1%.  Velocity  mea¬ 
surements  were  conducted  in  coincidence  mode  with 
20,000  random  samples  per  location.  Doppler  bursts  for 
the  two  channels  were  required  to  fall  within  a  50  ps 
coincidence  window  or  the  sample  was  rejected. 

In  this  study,  the  friction  velocity,  UT>  for  the  smooth 
surface  was  found  using  the  Clauser  chart  method  with 
log-law  constants  /c=0.41  and  B= 5.0.  For  the  rough  walls, 
UT  and  e  were  obtained  using  a  procedure  based  on  the 
modified  Clauser  chart  method  given  by  Perry  and  Li 
(1990).  For  the  rough  surfaces  the  top  plane  of  the 
roughness  elements  was  taken  to  be  y= 0.  The  wall  datum 
shift,  e,  is  the  distance  below  top  plane  of  the  roughness 
elements  that  U=0.  In  the  present  experiments,  e/Rt  gen¬ 
erally  ranged  between  0.25  and  0.35  for  the  roughness 
types  tested.  For  all  the  test  surfaces,  the  total  stress 
method  was  also  used  to  verify  Ux.  This  method  assumes  a 
constant  stress  region  equal  to  the  wall  shear  stress  exists 
in  the  inner  layer  of  the  boundary  layer.  If  the  viscous  and 
turbulent  stress  contributions  are  added  together,  an 
expression  for  UT  may  be  calculated  as  the  following 
evaluated  at  the  total  stress  plateau  in  the  inner  layer: 


In  all  cases  the  agreement  between  the  Clauser  chart 
and  total  stress  methods  was  well  within  their  uncertainty. 
The  roughness  function  for  the  rough  surfaces  was  ob¬ 
tained  as  the  difference  between  the  log-law  intercept  for 
the  smooth  and  the  rough  walls.  The  test  conditions  for  the 
velocity  profile  method  are  given  in  Table  2. 

2.2 

Rotating  disk  method 

The  rotating  disk  experiments  were  conducted  in  the  disk 
drag  facility  at  the  Naval  Surface  Warfare  Center, 
Carderock  Division.  A  schematic  of  the  facility  is  shown  in 
Fig.  2.  The  torque  measurements  were  made  on  22.9  cm 
diameter,  5  mm  thick  aluminum  disks,  coated  with  the 
specified  roughness,  rotated  in  the  enclosed  acrylic  tank. 
The  tank  is  25  cm  in  height  and  33  cm  in  diameter.  The 
angular  velocity  of  the  disks  was  varied  between  700  and 


Table  2.  Experimental  test 

conditions  for  velocity  profile  Specimen 

method 


Ue  Ree  CfXlO3  CfXlO3  <5  d*  H  Tl 

(ms-1)  Clauser  Total  stress  (mm)  (mm) 


Smooth 

0.94 

2,950 

3.44 

3.32 

28 

3.8 

1.30 

0.36 

2.60 

7,020 

2.99 

3.04 

26 

3.2 

1.27 

0.30 

2.99 

8,080 

2.92 

2.82 

27 

3.2 

1.26 

0.30 

3.58 

9,680 

2.82 

2.77 

26 

3.2 

1.26 

0.30 

60-Grit  sandpaper 

0.93 

3,720 

4.82 

4.55 

33 

5.1 

1.38 

0.40 

2.53 

10,600 

5.04 

5.29 

33 

5.5 

1.42 

0.38 

3.12 

13,800 

4.87 

5.09 

33 

5.9 

1.44 

0.44 

3.58 

16,400 

4.84 

5.13 

34 

6.1 

1.43 

0.45 

220-Grit  sandpaper 

0.95 

3,420 

3.52 

3.66 

33 

4.7 

1.36 

0.55 

2.60 

8,930 

3.79 

3.90 

29 

4.3 

1.34 

0.38 

3.07 

11,000 

3.89 

3.77 

30 

4.5 

1.36 

0.41 

3.63 

12,900 

3.85 

3.69 

30 

4.6 

1.36 

0.44 

Epoxy 

0.93 

3,170 

3.40 
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1,500  rpm  (ReR=9.0xl05-2.0xl06).  This  was  measured 
using  a  Scientific  Atlanta/Spectral  Dynamics  model  13134/ 
13135  GPT  stroboscopic  tachometer  (Spectral  Dynamics, 
San  Jose,  Calif.,  USA).  The  combined  precision  and  bias 
error  of  the  device  is  <1%  of  the  reading.  A  Baldor  model 
3445  DC  motor  (Baldor  Electric  Co.,  Fort  Smith,  Ariz., 
USA)  rated  at  1  hp  powers  the  rotating  disk  facility. 
Feedback  control  to  the  Baldor  model  BD154  motor  con¬ 
troller  was  used  to  hold  the  angular  velocity  constant 
based  on  input  from  the  tachometer.  The  torque  was 
measured  by  means  of  a  Lebow  model  1104-50  slip  ring- 
type  torque  sensor  (Lebow  Products,  Troy,  Mich.,  USA). 
This  sensor  has  a  range  of  0-5.7  N  m  with  a  combined  bias 
uncertainty  of  ±0.05%  of  full  scale.  The  working  fluid  in 
the  experiments  was  fresh  water,  and  the  temperature  was 
monitored  to  within  ±0.05°C  during  the  course  of  the 
experiments  using  a  thermocouple  logging  to  the  A/D 
system.  Data,  which  included  temperature,  angular  veloc¬ 
ity,  and  torque  were  gathered  at  a  sampling  rate  of 
1,000  Hz  and  were  digitized  using  a  16-bit  A/D  converter. 

Before  data  were  collected  on  the  surfaces  of  interest, 
the  tare  torque  values  caused  by  frictional  resistance  in  the 
bearings  and  hydrodynamic  torque  on  the  shaft  itself  were 
measured.  This  was  accomplished  by  rotating  the  shaft, 
without  the  disk,  in  the  tank  of  water.  This  was  done  by 
setting  the  angular  velocity  to  700  rpm  and  holding  it  there 


for  120  s.  Since  some  time  was  necessary  for  the  tank 
conditions  to  reach  an  equilibrium  state,  only  data  from 
the  final  30  s  were  used.  This  procedure  was  then  repeated 
for  900,  1,100,  1,300,  and  1,500  rpm.  Replicate  measure¬ 
ments  at  all  of  the  angular  velocities  were  made.  The  test 
disks  were  then  attached  to  the  shaft,  and  a  similar  process 
was  repeated  for  the  four  test  surfaces. 

It  should  be  noted  that  spinning  the  disk  in  the 
confined  tank  sets  up  a  mean  tangential  flow  in  the  fluid. 
The  result  is  that  the  actual  angular  velocity  of  the  disk 
relative  to  the  fluid  is  somewhat  less  than  the  angular 
velocity  of  the  disk  itself.  Granville  (1982)  states  that  this 
can  be  corrected  for  by  comparing  Cm  values  measured 
for  a  given  surface  in  the  tank  to  those  measured  for  the 
same  surface  in  a  very  large  body  of  water  (i.e.,  uncon¬ 
strained  fluid)  at  the  same  value  ReR(Cm)1/2.  He  terms 
this  the  swirl  factor  for  the  test  facility,  which  is  defined 
as  follows: 


The  test  facility  used  in  this  study  was  previously  used 
by  Loeb  et  al.  (1984),  who  found  0=0.854  over  a  wide 
range  of  ReR(Cm)1/2.  This  value  of  0  was,  therefore,  used  in 
the  present  study. 

The  roughness  functions  for  the  rotating  disk  data  were 
calculated  using  the  similarity  law  analysis  of  Granville 
(1982).  This  procedure  involves  comparing  the  Cm  values 
of  smooth  and  rough  disks  at  the  same  value  of  ReR(Cm)1/2. 
The  resulting  equations  for  k+  and  A  U1"  are  given  in 
Eqs.  (5)  and  (6),  respectively. 


Details  of  the  development  of  these  equations  are  given 
in  Granville  (1982).  It  should  be  noted  that  the  solutions  to 


both  of  the  equations  are  iterative.  Initially,  AU+'  is  taken 
to  be  zero  and  the  equations  are  solved.  Based  on  this 
solution,  AU+'  is  calculated  and  the  procedure  repeated 
until  the  solution  converges. 

2.3 

Towed  plate  method 

The  towed  plate  experiments  were  conducted  in  the  115-m 
long  towing  tank  facility  at  the  United  States  Naval 
Academy  Hydromechanics  Laboratory  using  an  experi¬ 
mental  method  similar  to  that  used  by  Schultz  (2002).  The 
width  and  depth  of  the  tank  are  7.9  m  and  4.9  m, 
respectively.  The  towing  carriage  has  a  velocity  range  of 
0-7.6  m/s.  In  the  present  investigation,  the  towing  velocity 
was  varied  between  2.0  m/s  and  3.8  m/s  (ReL=2.8xl06- 
5.5xl06).  The  velocity  of  the  towing  carriage  was  measured 
and  controlled  using  an  encoder  on  the  rails  that  produce 
4,000  pulses/m.  Using  this  system,  the  precision  uncer¬ 
tainty  in  the  mean  velocity  measurement  was  <0.02%  over 
the  entire  velocity  range  tested.  The  working  fluid  in  the 
experiments  was  fresh  water,  and  the  temperature  was 
monitored  to  within  ±0.05°C  during  the  course  of  the 
experiments  using  a  thermocouple  with  digital  readout. 

Figure  3  shows  a  schematic  of  the  test  fixture  and  plate. 
The  flat  test  plate  was  fabricated  from  304  stainless  steel 
sheet  stock  and  measured  1.52  m  in  length,  0.76  m  in 
width,  and  3.2  mm  in  thickness.  Both  the  leading  and 
trailing  edges  were  filleted  to  a  radius  of  1.6  mm.  No 
tripping  device  was  used  to  stimulate  transition.  The 
overall  drag  of  the  plate  was  measured  using  a  model 
HI-M-2,  modular  variable-reluctance  displacement  force 
transducer  manufactured  by  Hydronautics  Inc.  An  iden¬ 
tical  force  transducer,  rotated  90°  to  the  drag  gage,  was 
included  in  the  test  rig  to  measure  the  side  force  on  the 
plate.  The  purpose  of  the  side  force  gage  was  to  ensure 
precise  alignment  of  the  plate.  This  was  accomplished  by 
repeatedly  towing  the  plate  at  a  constant  velocity  and 
adjusting  the  yaw  angle  of  the  test  fixture  to  minimize  the 
side  force.  Once  this  was  done,  no  further  adjustments 
were  made  to  the  alignment  over  the  course  of  the 
experiments.  The  side  force  was  monitored  throughout  to 


TOWING  CARRIAGE  DIRECTION 


Fig.  3.  Schematic  of  the  flat  plate  test  fixture  for  the  towed  plate 
method 


confirm  that  the  plate  alignment  did  not  vary  between  test 
surfaces.  Both  of  the  force  transducers  used  in  the  exper¬ 
iments  had  load  ranges  of  0-89  N.  The  combined  bias 
uncertainty  of  the  gages  is  ±0.25%  of  full  scale.  Data  were 
gathered  at  a  sampling  rate  of  100  Hz  and  were  digitized 
using  a  16-bit  A/D  converter.  The  length  of  the  towing  tank 
dictated  the  sampling  duration.  This  ranged  from  ~30  s  of 
data  per  test  run  at  the  lowest  Reynolds  number  to  ~11  s 
of  at  the  highest  Reynolds  number.  The  overall  drag  was 
first  measured  with  590  mm  of  the  plate  submerged.  This 
was  repeated  with  25  mm  of  the  plate  submerged  in  order 
to  find  the  wavemaking  resistance  tare.  The  difference 
between  the  two  was  taken  to  be  the  frictional  resistance 
on  the  two  565  mm  wide  by  1.52  m  long  faces  of  the  plate. 

The  roughness  functions  for  the  towed  plate  data  were 
calculated  using  the  similarity  law  analysis  of  Granville 
(1987).  This  procedure  involves  comparing  the  CF  values 
of  smooth  and  rough  plates  at  the  same  value  of  ReLCF. 
The  resulting  equations  for  k+  and  A  U+  are  given  in 
Eqs.  (7)  and  (8),  respectively. 


Details  of  the  development  of  these  equations  are  given 
in  Granville  (1987).  Again,  the  solutions  to  both  of  the 
equations  are  iterative.  Initially,  A  U+'  is  taken  to  be  zero 
and  the  equations  are  solved.  Based  on  this  solution,  A  U+' 
is  calculated  and  the  procedure  repeated  until  the  solution 
converges. 

2.4 

Uncertainty  estimates 

Precision  uncertainty  estimates  for  all  of  the  measure¬ 
ments  in  this  investigation  were  made  through  repeat¬ 
ability  tests.  This  involved  conducting  repeated 
experiments  for  a  given  surface  condition  using  each  of  the 
testing  methods.  In  order  to  estimate  the  95%  precision 
confidence  limits  for  a  measured  parameter,  its  standard 
deviation  was  multiplied  by  the  two-tailed  t  value  for 
(N-l)  degrees  of  freedom,  as  given  by  Coleman  and  Steele 
(1995).  These  uncertainties  were  then  combined  with  the 
estimated  bias  uncertainties  using  standard  error  propa¬ 
gation  methods  (Moffat  1988)  to  find  the  overall  95% 
confidence  limits  for  the  calculated  parameters.  The 
resulting  overall  95%  confidence  limits  for  Cm  and  CF  were 
both  ±2%.  The  resulting  overall  95%  confidence  limits  for 
Cf  were  ±4%  for  the  smooth  wall  and  ±7%  for  the  rough 
walls  using  the  Clauser  chart  method.  The  Cf  results  using 
the  total  stress  method  (Eq.  3)  agreed  within  their  uncer¬ 
tainty  with  those  obtained  using  the  Clauser  chart  method. 
The  resulting  uncertainties  in  A U+  were  ±10%  or  0.2 
(whichever  was  greater),  ±9%  or  0.3  (whichever  was 
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Fig.  4.  Skin-friction  coefficient  vs.  momentum  thickness  Reynolds  Fig.  6.  Torque  coefficient  vs.  Reynolds  number  for  the  rotating 
number  for  the  velocity  profile  method.  (Overall  uncertainty  in  disk  method.  (Overall  uncertainty  in  Cm:  ±2%) 
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Fig.  5.  Mean  velocity  profiles  in  wall  coordinates  for  all  surfaces  at 
the  highest  Rex.  (Overall  uncertainty  in  CJ' :  smooth  wall,  ±3%; 
rough  wall,  ±5%) 


greater),  and  ±12%  or  0.2  (whichever  was  greater)  for  the 
velocity  profile,  rotating  disk,  and  towed  plate  methods, 
respectively. 

3 

Results  and  discussion 

The  skin-friction  results  for  the  velocity  profile  method  are 
shown  in  Fig.  4.  The  smooth  wall  results  of  Coles  (1962) 
are  presented  for  comparison.  Both  sandpaper  rough 
surfaces  exhibited  a  significant  increase  in  Cf  over  the 
entire  range  of  Ree.  For  example,  at  the  highest  Reynolds 
number,  Cf  was  86%  higher  for  the  60-grit  sandpaper  and 
was  42%  higher  for  the  220-grit  sandpaper  than  the 
smooth  wall  values.  The  Cf  values  for  epoxy  surface  are 
observed  to  rise  only  slightly  above  the  smooth  wall  values 
as  Reo  increases.  At  the  highest  Reynolds  number,  the  Cf 
for  the  epoxy  surface  is  8%  higher  than  the  smooth  wall 
value,  which  is  still  within  the  combined  uncertainty  of  the 
measurements. 

Figure  5  shows  the  mean  velocity  profiles  in  wall 
variables  for  all  of  the  test  surfaces  at  the  highest  Rex  tested 
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Fig.  7.  Frictional  resistance  coefficient  vs.  Reynolds  number  for 
the  towed  plate  method.  (Overall  uncertainty  in  CF:  ±2%) 


(~4.9xl06).  The  results  show  that  the  smooth  profile  fol¬ 
lows  the  classical  log-law  (Eq.  1)  well  in  the  overlap  region. 
The  three  rough  surfaces  also  display  a  linear  log  region 
that  is  shifted  by  AC1-  below  the  smooth  profile  (Eq.  2).  A 
trend  of  increasing  AU+  with  increasing  roughness  height 
is  observed  in  the  results. 

Figure  6  shows  the  results  of  the  rotating  disk  experi¬ 
ments.  Shown  for  comparison  is  the  empirical  relationship 
of  von  Karman  (1946)  for  the  torque  coefficient  of  an 
unconfined,  smooth  rotating  disk  given  as: 


(9) 


The  smooth  disk  Cm  results  show  good  agreement 
(within  ±1%)  with  Eq.  9  over  the  entire  range  of  ReR.  All  of 
the  rough  disks  exhibited  significant  increases  in  Cm.  The 
epoxy-coated  disk  had  an  increase  in  Cm  of  22-37% 
compared  to  the  smooth  disk  values.  The  increases  in  Cm 
for  the  220-grit  and  60-grit  sandpaper  surfaces  were 
88-116%  and  150-179%,  respectively.  The  reason  for  the 
larger  increases  observed  in  Cm  compared  to  Cf  on  the 
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Fig.  8.  Roughness  function  for  the  test  surfaces  obtained  using  the 
three  methods.  (Overall  uncertainty  in  Al/+:  velocity  profile 
method,  ±10%  or  0.2  (whichever  is  greater);  rotating  disk 
method,  +9%  or  0.3  (whichever  is  greater);  towed  plate  method, 
±12%  or  0.2  (whichever  is  greater)) 

same  surface  was  that  the  rotating  disk  experiments  were 
conducted  at  unit  Reynolds  numbers  of  up  to  five  times 
higher  than  the  velocity  profile  experiments.  This  leads  to 
an  increased  effect  for  a  roughness  of  a  given  height. 

Figure  7  shows  the  results  of  the  towed  plate  experi¬ 
ments.  The  Karman-Schoenherr  friction  line  (Schoenherr 
1932)  for  the  frictional  resistance  of  a  smooth  flat  plate  is 
also  presented.  It  is  given  as  follows: 

0-242  ,  , 

-j=  =  log  (ReLCv)  (10) 

The  smooth  plate  CF  values  agree  within  ±2%  with  the 
Karman-Schoenherr  friction  line.  The  CF  values  for  the 
220-grit  and  60-grit  sandpaper  surfaces  were  increased 
significantly  compared  to  the  smooth  plate  values  over  the 
entire  Rex  range.  The  increases  in  CF  for  the  220-grit  and 
60-grit  sandpaper  surfaces  were  17-30%  and  66-83%, 
respectively.  The  CF  values  for  the  epoxy-coated  plate 
showed  only  modest  increase  compared  to  the  smooth 
plate.  At  the  highest  Rex,  the  increase  was  6%,  which  is 
slightly  greater  than  the  combined  uncertainty  in  the 
measurements. 

The  roughness  functions  for  each  of  the  test  surfaces 
using  the  three  determination  methods  are  presented  in 
Fig.  8.  Shown  for  comparison  is  the  Nikuradse-type 
roughness  function  for  uniform  sand  given  by  Schlichting 
(1979).  Also  shown  is  the  Colebrook-type  roughness 
function  for  random  roughness  given  by  Grigson  (1992)  as 
follows: 

AU+  =  -  In  (1  +  k+)  (11) 

K 

The  chosen  roughness  length  scale  for  the  sandpaper 
surfaces  was  k=0.75Rt.  This  was  the  length  scale  suggested 
by  Schultz  and  Flack  (2003)  for  sandpaper  roughness.  The 
roughness  length  scale  used  for  the  epoxy  surfaces  was 
k=039Ra.  This  length  scale  was  offered  by  Schultz  and 
Flack  (2003)  as  an  appropriate  one  for  a  range  of  painted 
and  sanded  surfaces.  It  should  be  noted  that  the  choice  of 


k  for  a  given  roughness  has  no  effect  on  the  calculated  ALT' 
despite  its  apparent  dependence  on  k  through  the  ALL1"' 
term  in  Eqs.  6  and  8.  This  is  because  the  effect  of  changing 
k  on  ALL1"  is  to  simply  move  the  curve  along  the  abscissa 
without  changing  its  slope. 

The  roughness  functions  obtained  using  the  velocity 
profile  and  towed  plate  methods  show  good  agreement  for 
all  three  rough  surfaces.  The  epoxy  surface  follows  a 
Colebrook-type  roughness  function  using  the  roughness 
length  scale  used  here.  The  sandpaper  roughness  results 
using  these  methods  agree  within  their  uncertainty  with  a 
Nikuradse-type  roughness  function.  Since  the  rotating 
disk  results  were  obtained  at  much  higher  unit  Reynolds 
numbers  than  the  other  two  methods,  direct  comparison 
of  the  results  is  not  possible.  However,  using  the  rotating 
disk  method,  the  epoxy  surface  again  has  ALT'  values  that 
agree  within  their  uncertainty  with  a  Colebrook-type 
roughness  function.  The  sandpaper  surfaces  have  AlT" 
values  that  agree  within  uncertainty  with  a  Nikuradse-type 
roughness  function.  Furthermore,  the  experimental  results 
for  both  sandpaper  surfaces  have  slopes  (ALT')  of  ~1  Ik, 
the  slope  the  roughness  function  is  expected  to  exhibit  in 
the  fully  rough  flow  regime. 

4 

Conclusion 

Three  experimental  methods  have  been  used  to  determine 
the  roughness  function  of  several  surfaces.  The  results  for 
the  velocity  profile  and  towed  plate  methods  show  good 
agreement  for  all  surfaces,  with  the  epoxy  surface  follow¬ 
ing  a  Colebrook-type  roughness  function  and  the  sand¬ 
paper  surfaces  following  a  Nikuradse-type  roughness 
function.  Although  the  tests  using  the  rotating  disk 
method  were  carried  out  at  much  higher  unit  Reynolds 
numbers,  the  results  for  the  sandpaper  rough  surfaces 
agree  within  their  uncertainty  with  a  Nikuradse-type 
roughness  function  in  the  fully  rough  regime,  while  the 
results  for  the  epoxy  surface  agree  with  a  Colebrook-type 
roughness  function.  The  present  results  illustrate  the 
utility  of  indirect  methods  in  roughness  function  deter¬ 
mination  and  also  indicate  that  the  rotating  disk  method 
may  provide  a  relatively  economical  means  for  determin¬ 
ing  the  effect  of  surface  roughness  at  high  Reynolds 
number. 
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